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• Premise of the study: Three species of the mangrove tree genus Rhizophora are found in the New World and along the west 
coast of Africa. Of these, R. mangle is the most abundant and has a complex interbreeding relationship with the sympatric 
R. racemosa and R. harrisonii. The development of additional microsatellite markers would permit paternity analyses and 
investigation of the hybrid origin of these species. 

• Methods and Results: Using an enriched library method, via hybridization with biotinylated oligonucleotides complementary 
to repetitive poly AG/TC, primers for 1 1 microsatellite markers of R. mangle were developed and characterized in populations 
in Para and Sao Paulo (Brazil) and Florida (USA). Ten of these markers were transferable to R. racemosa and R. harrisonii. 

• Conclusions: The microsatellite markers presented here will be useful in studies of contemporary and historical gene flow be- 
tween American and West African Rhizophora species. 
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In the New World and along the west coast of Africa, three 
species of the mangrove tree genus Rhizophora L. (Rhizopho- 
raceae) have been identified. Of these, R. mangle L. (red man- 
grove) is the most widely distributed (Tomlinson, 1986). In 
terms of interbreeding, R. mangle has a complex relationship 
with the sympatric R. racemosa G. Mey. and R. harrisonii Leechm. 
(Ceron-Souza et al., 2010), which could potentially be resolved 
by the use of microsatellite markers. Two studies have pre- 
viously reported the isolation of microsatellite markers for 
R. mangle. In Rosero-Galindo et al. (2002), 10 polymorphic 
loci with between two and seven alleles per locus were reported 
(26 alleles total). In Takayama et al. (2008), 14 simple sequence 
repeat (SSR) loci were reported and 1 1 were shown to be poly- 
morphic, with a total of 32 and 27 alleles in populations from 
the Pacific and Atlantic coasts of Costa Rica, respectively. The 
allelic diversity reported for previously isolated microsatellite 
markers for this species is therefore somewhat low, and this is 
also the case in subsequent studies that made use of these mark- 
ers. Arbelaez-Cortes et al. (2007) report a total of 17 alleles 
for three markers in five studied populations, Pil et al. (2011) 
report 22 alleles for eight markers in 10 populations, and 
Sandoval-Castro et al. (2012) report 19 alleles for six markers 
in 10 populations. For this reason, we isolated and characterized 
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new microsatellite markers for R. mangle, which were then sub- 
sequently tested for transferability to R. racemosa and R. har- 
risonii. These markers present higher allelic diversity than has 
previously been reported. In conjunction with pre-existing 
markers, these new markers will permit mating systems analysis, 
paternity tests, and future studies attempting to resolve the 
potential hybrid origin of Rhizophora species. 



METHODS AND RESULTS 

Genomic DNA was isolated from leaf tissue of one adult individual of 
R. mangle sampled in the Caete estuary, Braganca, Para, northern Brazil 
(0°52 , 22 ,, S, 46°39 , 04"W), and used in the development of the enriched library. 
Herbarium material was deposited at the herbarium of the Institute for Coastal 
Research in Braganca (HBRA; accession no. 821). Polymorphism was evalu- 
ated using the 24 adult trees of R. mangle from Para, Brazil (0°50'27"S, 
46°38'15"W; HBRA 706); Sao Paulo, Brazil (24°00 , 21 ,, S, 46°17 , 59"W); and 
Florida, USA (25°08'09"N, 80°57'29"W; HBRA 1291). The Sao Paulo and 
Florida sites are near the southernmost and northernmost distributional limits 
of R. mangle, respectively. For the transferability test, we genotyped seven in- 
dividuals of R. harrisonii and five individuals of R. racemosa. 

An enriched library was constructed according to Vinson et al. (2005). Total 
DNA (50,000 ng) was digested with Sau31 and separated on a 2% agarose gel. 
Fragments of 300-800 bp were recovered with purification of the gel using 
QIAquick Gel Extraction kit (QIAGEN, Sao Paulo, Brazil). DNA fragments 
were ligated to adapters, hybridized to biotinylated (AG) 13 and (TC) 13 , and the 
oligofragments separated using streptavidin magnetic beads. Selected frag- 
ments were ligated into a pGEM-T Easy Vector (Promega Corporation, Madison, 
Wisconsin, USA) and transformed in competent E. coli XL 1 -Blue cells that 
were grown overnight on lx Luria-Bertani agar plates containing ampicillin, 
X-gal, and IPTG. Transformants (white colonies) were diluted in water, and the 
fragment inserted into the plasmid was amplified using the Ml 3 universal prim- 
ers (M13F[-20] and M13R[-40]) by PCR. PCR products were sequenced on a 
MegaBACE 1000 instrument (GE Healthcare, Belo Horizonte, Minas Gerais, 
Brazil) using dye terminator fluorescent chemistry (Sanger method) with the 
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Table 1. Characteristics of 1 1 microsatellite loci developed in Rhizophora mangle. 



Loci 




Primer sequences (5'-3') 


Repeat motif 


Allele size (bp) 


T a (°Q 


GenBank accession no. 


RmBral8 


F 
R 


TCAAGGACAGGTCAACAGCA 
TGAATTAAGACGTCAAATCATCG 


(TC) 18 


174 


54 


HQ850264 


RmBral9 


F 
R 


GAGGCAGAGTCAGGTCAGAA 
CACTGGTCCACTGACAGCAA 


(CT) 14 


162 


54 


HQ850265 


RmBra20 


F 
R 


TCAGCACAATACATCAGGACAA 
GGCGCTATCCATCCTGAGT 


(AG) 16 


141-175 


50 


HQ850255 


RmBra25 


F 
R 


TCCATTCTGCCAAATTTGATT 
CGAATTCACCGATCACAATG 


(AG) 13 


172-226 


50 


HQ850256 


RmBra27 


F 
R 


TCTTTTGTCGAAGCCTCTCC 
ATCAAGTTGGCGTTGGAAAA 


(CT) 3 (TC) 6 (CT) 4 (TC) 3 


140-184 


50 


HQ850257 


RmBra45 


F 
R 


GAAAATGCAAGAGGGCTGAC 
CATGGGTCTCTCTCCTCTGC 


(AG) 17 (AGG) 3 


128-164 


54 


HQ850258 


RmBra50 


F 

R 


ATCGTGGAAGAACGGGGTTT 
TCAAGAAGTCCAGGGTGCTT 


(AG) 18 


149-215 


50 


HQ850259 


RmBra59 


F 

R 


GTGAACGGTCTGGACTGGAG 
TCACCGATCCCCTAGAACTG 


(AG) 24 


145-209 


54 


HQ850260 


RmBra64 


F 

R 


CAAGGCTAATCCGAAAATCG 
ACTCTGTCCCCTGCCTCAC 


(AG) 15 


148 


54 


HQ850266 


RmBra65 


F 

R 


CCATAGACAATACAGGATACCCAGA 
CAGCATGAACATCACCTTGG 


(AG) 12 


234 


52 


HQ850263 


RmBra66 


F 

R 


TCAGCACAATACATCAGGACAA 
GCGCTATCCATCCTGAGTTT 


(AG) 19 


138-206 


54 


HQ850261 



Note: T a = annealing temperature. 



M13F(-20) primer. The sequences were subsequently edited manually using 
Sequencher version 4.1.4 software (Gene Codes Corporation, Ann Arbor, 
Michigan, USA). Primers were designed in the regions flanking AG/CT repeats 
using Primer3 OUTPUT (Rozen and Skaletsky, 2000) with the following condi- 
tions: final amplicon between 100 and 350 bp, average annealing temperature 
of 55°C, primer length between 18 and 22 bp, and GC content between 40% 
and 60%. We sequenced 66 clones, and 23 pairs of primers were synthesized. 
The preliminary polymorphism test was conducted with 12 individuals from the 
Caete estuary (Braganca). PCR products were visually compared with the 
10-bp ladder in a 10% polyacrylamide gel stained with ethidium bromide. 

From the 23 primers, 1 1 primers produced PCR products, seven loci were 
polymorphic, and four loci were monomorphic (RmBral8, RmBral9, RmBra64, 
and RmBra65; Table 1). Forward primers of the seven loci were fluorescently 
labeled with 6-FAM (RmBra20, RmBra25, and RmBra59; MWG-Biotech, 
Ebersberg, Germany) or HEX (RmBra27, RmBra45, RmBra50, and RmBra66; 
MWG-Biotech). Loci were amplified using PCR in a total volume of 13 uL 
containing: DNA (5 ng), lx PCR reaction buffer (10 mM Tris-HCl [pH 8.3], 
50 mM KC1, 1.5 mM MgCl 2 ), forward primer and reverse primer (0.15 jiM), 
MgCl 2 (1.5 mM), bovine serum albumin (BSA; New England Biolabs, Hitchin, 
Hertfordshire, United Kingdom; 0.75 mg-mL" 1 ), dNTP (0.5 mM), and 1.3 U Taq 
polymerase (Invitrogen, Life Technologies, Carlsbad, California, USA). PCR 
conditions were: denaturation at 96°C for 2 min; 30 cycles of denaturation at 
94°C for 15 s, annealing at 50-54°C (Table 1) for 15 s, and extension at 72°C 
for 15 s; and a final extension at 72°C for 5 min. Transfer of the primers of 
R. mangle to R. racemosa and R. harrisonii was tested under the same con- 
ditions as for R. mangle (described above). Fragments were genotyped on a 
96-capillary sequencer MegaBACE 1000 (GE Healthcare). Fragments were 
analyzed using MegaBACE Fragment Profiler version 1 .2 (Amersham Biosci- 
ence, Belo Horizonte, Minas Gerais, Brazil) using MegaBACE ET 550-R Rox 
Size Standard (Amersham Bioscience). Number of alleles per locus, observed 
and expected heterozygosities, coefficient of fixation index, and Hardy-Weinberg 
and linkage disequilibrium were calculated using FSTAT (Goudet, 2002). Probabil- 
ities of paternity exclusion were estimated using CERVUS (Kalinowski et al., 
2007). Frequency of null alleles was estimated using MICRO-CHECKER version 
2.2.3 (van Oosterhout et al., 2004). 

Results show a total of 90 alleles in the three studied populations (49 in 
Braganca, 39 in Sao Paulo, and 42 in Florida), and the number of private alleles 
varied between one to eight alleles per locus (Table 2). Observed heterozygosity 
was lower than values of heterozygosity under expected Hardy-Weinberg equi- 
librium (HWE) for all populations tested, with higher fixation index values for 
the Sao Paulo and Florida populations (Table 2). The test of pairwise linkage 
after the Bonferroni correction showed values higher than the significant value 
of 0.005, indicating that these loci are not linked and can be used as indepen- 
dent markers. Only one locus, RmBra27, showed significant null alleles using 



the Oosterhout index (0.26). Fisher's exact test revealed deviations from HWE 
after Bonferroni correction at two loci for all populations tested, with an excess 
of homozygotes for the loci RmBra25, RmBra27, and RmBra50 (P = 0.019). 
The combined exclusion power of all seven loci was 0.998, 0.998, and 0.982 for 
first parent total exclusion probabilities for the Braganca, Sao Paulo, and Florida 
populations, respectively. Transfer of the primers was successful, with six poly- 
morphic markers identified for R. mangle also polymorphic in R. racemosa and 
R. harrisonii (except RmBra25), while the four monomorphic markers identi- 
fied for R. mangle were only tested for amplification. 



CONCLUSIONS 

The high number of private alleles in the studied populations 
indicates that these loci are useful for studies of genetic divergence, 
genetic structure, and phylogeography. In addition, the parental 
exclusion probability values are adequate for progeny analysis 
and paternity tests for further identification of reproduction bar- 
riers and identification of hybridization between species of all 
American and West African Rhizophora species. 
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